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METHOD AND APPARATUS FOR STORING BITPLANES OF 
COEFFICIENTS IN A REDUCED SIZE MEMORY 

FIELD OF THE INVENTION 

[0001] The present invention relates to the field of compression and 
decompression; more particularly, the present invention relates to storing 
bitplanes of coefficients in a reduced size memory. 

BACKGROUND OF THE INVENTION 

[0002] Data compression is an extremely useful tool for storing and 
transmitting large amotints of data. For example, the time required to 
transmit an image, such as a facsimile transmission of a document, is 
reduced drastically when compression is used to decrease the number of bits 
required to represent the image. 

[0003] Many different data compression techniques exist in the prior art. 
Compression techniques can be divided into two broad categories, lossy 
coding and lossless coding. Lossy coding involves coding that results in the 
loss of information, such that there is no guarantee of perfect reconstruction 
of the original data. The goal of lossy compression is that changes to the 
original data are done in such a way that they are not objectionable or 
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detectable. In lossless compression, all the information is retained and the 
data is compressed in a manner which allows for perfect reconstruction. 
[00041 In lossless compression, input symbols or intensity data are converted 
to output codewords. The input may include image, audio, one-dimensional 
(e.g., data changing spatially), two-dimensional (e.g., data changing in two 
spatial directions), or multi-dimensional/multi-spectral data. If the 
compression is successful, the codewords are represented in fewer bits than 
the number of bits in the "normal" representation of the input symbols (or 
intensity data). Lossless coding methods include dictionary methods of 
coding (e.g., Lempel-Ziv), rim length encoding, enumerative coding and 
entropy coding. In lossless image compression, compression is based on 
predictions or contexts, plus coding. The JBIG standard for facsimile 
compression and DPCM (differential pulse code modulation - an option in 
the JPEG standard) for continuous-tone images are examples of lossless 
compression for images. In lossy compression, input symbols or intensity 
data are quantized prior to conversion to output codewords. Quantization 
is intended to preserve relevant characteristics of the data while eliminating 
xmimportant characteristics. Prior to quantization, lossy compression 
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system often use a transform to provide energy compaction. JPEG is an 
example of a lossy coding method for image data. 

[0005] Reversible transforms (wavelet, component) may be used for both 
lossy and lossless compression. Irreversible transforms (wavelet, 
component, discrete cosine) may be used only for lossy. 
[0006] The new JPEG 2000 decoding standard utilizes transforms and 
provides a new coding scheme and codestream definition for images. 
Although JPEG 2000 is a decoding standard, and thus defines what a 
decoder must do, this definition restricts an encoder especially for lossless 
compression. Under the JPEG 2000 Standard, each image may be divided 
into rectangular tiles. If there is more than one tile, the tiling of the image 
creates tile-components. An image may have multiple components. For 
example, a color image might have red, green and blue components. Tile- 
components can be extracted or decoded independently of each other. 
[0007] After tiling of an image, the tile-components are decomposed into one 
or more different decomposition levels using a wavelet transformation. 
These decomposition levels contain a number of subbands populated with 
coefficients that describe the horizontal and vertical spatial frequency 
characteristics of the original tile-components. The coefficients provide 
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frequency information about a local area, rather than across the entire image. 
That is, a small number of coefficients completely describe a single sample, 
A decomposition level is related to the next decomposition level by a spatial 
factor of two, such that each successive decomposition level of the subbands 
has approximately half the horizontal resolution and half the vertical 
resolution of the previous decomposition level. 

[0008] Although there are as many coefficients as there are samples, the 
information content tends to be concentrated in just a few coefficients. 
Through quantization, the numerical precision of a number of coefficients 
may be reduced with a disproportionately low introduction of distortion 
(quantization noise). Additional processing by an entropy coder reduces the 
number of bits required to represent these quantized coefficients, sometimes 
significantly compared to the original image. 

[0009] The individual subbands of a tile-component are further divided into 
code-blocks. These code-blocks can be grouped into precincts. These 
rectangular arrays of coefficients can be extracted independently. The 
individual bit-planes of the coefficients in a code-block are entropy coded 
with three coding passes. Each of these coding passes collects contextual 
information about the bit-plane compressed image data. 
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[0010] The bit stream compressed image data created from these coding 
passes is grouped in layers. Layers are arbitrary groupings of successive 
coding passes from code-blocks. Although there is great flexibility in 
layering, the premise is that each successive layer contributes to a higher 
quality image. Code-blocks of subband coefficients at each resolution level 
are partitioned into rectangular areas called precincts* 
[0011] Packets are a fundamental unit of the compressed codestream. A 
packet contains compressed image data from one layer of a precinct of one 
resolution level of one tile-component. These packets are placed in a defined 
order in the codestream. 

[0012] The codestream relating to a tile, organized in packets, are arranged 
in one, or more, tile-parts, A tile-part header, comprised of a series of 
markers and marker segments, or tags, contains information about the 
various mechanisms and coding styles that are needed to locate, extract, 
decode, and reconstruct every tile-component. At the begiiming of the 
entire codestream is a main header, comprised of markers and marker 
segments, that offers similar information as well as information about the 
original image. 
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[0013] The codestream is optionally wrapped in a file format that allows 
applications to interpret the meaning of, and other information about, the 
image. The file format may contain data besides the codestream. 
[0014] The decoding of a JPEG 2000 codestream is performed by reversing 
the order of the encoding steps. Figure 1 is a block diagram of the JPEG 2000 
standard decoding scheme that operates on a compressed image data 
codestream. Referring to Figure 1, a bitstream initially is received by data 
ordering block 101 that regroups layers and subband coefficients. 
Arithmetic coder 102 uses contextual information from previously coded 
coefficients provided by the bit modeling block 103 about the bit-plane 
compressed image data, and its internal state, to decode a compressed bit 
stream. 

[0015] Next, the codestream is quantized by quantization block 104, which 
may be quantizing based on a region of interest (ROI) as indicated by ROI 
block 105. After quantization, an inverse wavelet/spatial transform is 
applied to the coefficients via transform block 107, followed by DC level 
shifting and optional component transform block 108. This results in 
generation of a reconstructed image. 
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SUMMARY OF THE INVENTION 

[0016] A method and apparatus for storing bitplanes of coefficients in a 
reduced size memory is disclosed. In one embodiment, a method for 
encoding coefficients comprises applying one or more wavelet transforms to 
generate multiple coefficients and storing N bitplanes of each of the multiple 
coefficients in a first memory having multiple rows. Each of the multiple 
coefficients has a bit length greater than N bits, and each storage location in 
each of the rows stores a bit of either a more significant bitplane of the 
coefficients or a less significant bitplane of the coefficients. The storing of N 
bitplanes includes storing an indication for each row of the rows to indicate 
a location in each row at which bits of more significant bitplanes have 
started to be stored and after which bits of the multiple coefficients of the 
less significant bitplanes are no longer stored. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0017] The present invention will be understood more fully from the 
detailed description given below and from the accompanying drawings of 
various embodiments of the invention, which, however, should not be taken 
to limit the invention to the specific embodiments, but are for explanation 
and imderstanding only. 

[0018] Figure 1 is a block diagram of the JPEG 2000 standard decoding 
scheme. 

[0019] Figure 2A shows an example 8x8 code-block of coefficients with the 
sub-bit-plane pass identified for each coefficient and a label showing the 
order of processing for each coding pass. 

[0020] Figure 2B illustrates a memory for a variable length rxm and skip 
counts. 



[0021] Figures 3A-D illustrate neighborhood coefficients and memory 
organization for one embodiment of a context model. 
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[0022] Figure 4 shows one embodiment of a significance memory 
organization for random access of a 16x16 code-block. 

[0023] Figure 5 shows the memories and registers used in the significance 
propagation pass for random access. 

[0024] Figure 6 illustrates the significance state from memory and stored in 
registers for random access. 

[0025] Figure 7 is a block diagram of one embodiment of significance 
propagation pass logic. 

[0026] Figure 8 illustrates an example of performance of one embodiment of 
a context model on a 4x4 block. 

[0027] Figure 9 illustrates one embodiment of an organization of a 
significance memory for sequential accesses of a 16x16 code-block. 
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[0028] Figure 10 illustrates one embodiment of the memories and registers 
used for the significance propagation pass. 

[0029] Figure 11 shows how memory and registers may be used to provide 
the proper region for context model operation. 

[0030] Figure 12 is a block diagram of one embodiment of pass 
determination logic. 

[0031] Figure 13 is a blocl< diagram of one embodiment of double context 
generation logic. 

[0032] Figure 14A is a bloclc diagram of one embodiment of an "early 
context'' MQ-decoder. 

[0033] Figure 14B is one embodiment of a typical decoding implementation. 

[0034] Figure 15 is a block diagram of one embodiment of a "late context" 
MQ-decoder. 
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[0035] Figure 16A illustrates how a comparison of probability class indices 
may operate. 

[0036] Figure 16B is a block diagram of a multiplexor that determines the 
MPS or LPS for each Qe_Value. 

[0037] Figure 17 illustrates the intervals for multiple MPS decoding. 

[0038] Figure 18 illustrates one embodiment of an assignment of code blocks 
in parallel for 4:4:4 data. 

[0039] Figure 19 illustrates one embodiment of an assignment of code blocks 
in parallel for 4:2:2 data. 

[0040] Figure 20 illustrates an alternative embodiment of an assignment of 
code blocks in parallel for 4:1:1 data. 
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[0041] Figure 21 is a diagram of memory for one embodiment of a coder that 
includes multiple MQ coders, each having an associated context model 

[0042] Figure 22A shows how to use a limited number of bitplanes of 
memory for each coefficient during encoding. 

[0043] Figure 22B illustrates a bloc]< diagram of one embodiment of control 
logic to control access to the memory. 

[0044] Figure 23 illustrates using a small amoimt of memory for variable 
length (VL) code information instead of coxmters. 

[0045] Figure 24 is a block diagram of one embodiment of a forward 
transform. 

[0046] Figures 25A and 26B are block diagrams of embodiments of the low 
pass filter. 

[0047] Figures 26A and 25B are block diagrams of embodiments of a high 
pass filter. 
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[0048] Figure 27 is a block diagram of an alternative embodiment of a 
forward transform. 

[0049] Figure 28 is a block diagram of one embodiment of an inverse 
transform. 

[0050] Figure 29 is a block diagram of one embodiment of an 
encoder / decoder. 

[0051] Figure 30 illustrates an example of a 16 bit word having both run 
counts and skip counts. 

[0052] Figure 31 illustrates an exemplary 8x8 region of significance state bits 
for determining the coding pass. 

[0053] Figure 32 is a diagram of one embodiment of determine pass logic. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 
[0054] Techniques for performing coding are described. These techniques 
may be used for implementing JPEG 2000 or for manipulating or adding to 
the feature set. That is, the JPEG2000 standard, Information Technology - 
JPEG 2000 image coding system: Core Coding System, ITU-T Rec. 
T.800 1 ISO/IEC FDIS 15444-1: 2000 JPEG Image Coding System, 
incorporated herein by reference, leaves many choices to implementers. It is 
a purpose of the techniques described herein to use choices in JPEG 2000 to 
make high speed, low cost, low memory and/or feature rich 
implementations in software, hardware, and/ or firmware. 
[0055] In the following description, numerous details are set forth to provide 
a thorough xmderstanding of the present invention. It will be apparent, 
however, to one skilled in the art, that the present invention may be 
practiced without these specific details. In other instances, well-known 
structures and devices are shown in block diagram form, rather than in 
detail, in order to avoid obscuring the present invention. Furthermore, any 
block, logic or functionality not described in detail may be implemented 
with well-known implementations or may be easily implemented by those 
skilled in the art using well-known hardware, software and/or firmware. 
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Note that sometimes techniques and implementations are described using 
psuedo code. This is not meant to imply that the techniques are 
implemented solely in software; instead, such a description is often chosen 
to easily describe the functionality of terms one skilled in the art would 
easily imderstand* 

[0056] Some portions of the detailed descriptions which follow are presented 
in terms of algorithms and symbolic representations of operations on data 
bits within a computer memory. These algorithmic descriptions and 
representations are the means used by those skilled in the data processing 
arts to most effectively convey the substance of their work to others skilled 
in the art. An algorithm is here, and generally, conceived to be a self- 
consistent sequence of steps leading to a desired result. The steps are those 
requiring physical manipulations of physical quantities. Usually, though 
not necessarily, these quantities take the form of electrical or magnetic 
signals capable of being stored, transferred, combined, compared, and 
otherwise manipulated. It has proven convenient at times, principally for 
reasons of common usage, to refer to these signals as bits, values, elements, 
symbols, characters, terms, numbers, or the like. 
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[0057] It should be borne in mind, however, that all of these and similar 
terms are to be associated with the appropriate physical quantities and are 
merely convenient labels applied to these quantities. Unless specifically 
stated otherwise as apparent from the following discussion, it is appreciated 
that throughout the description, discussions utilizing terms such as 
"processing" or "computing" or "calculating" or "determining" or 
"displaying" or the like, refer to the action and processes of a computer 
system, or similar electronic computing device, that manipulates and 
transforms data represented as physical (electronic) quantities within the 
computer system's registers and memories into other data similarly 
represented as physical quantities within the computer system memories or 
registers or other such information storage, transmission or display devices. 
[0058] The present invention also relates to apparatus for performing the 
operations herein. This apparatus may be specially constructed for the 
required purposes, or it may comprise a general purpose computer 
selectively activated or reconfigured by a computer program stored in the 
computer. Such a computer program may be stored in a computer readable 
storage medium, such as, but is not limited to, any type of disk including 
floppy disks, optical disks, CD-ROMs, and magnetic-optical disks, read-only 
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memories (ROMs), random access memories (RAMs), EPROMs, EEPROMs, 
magnetic or optical cards, or any type of media suitable for storing electronic 
instructions, and each coupled to a computer system bus. 
[0059] The algorithms and displays presented herein are not inherently 
related to any particular computer or other apparatus. Various general 
purpose systems may be used with programs in accordance with the 
teachings herein, or it may prove convenient to construct more specialized 
apparatus to perform the required method steps. The required structure for 
a variety of these systems will appear from the description below. In 
addition, the present invention is not described with reference to any 
particular programming language. It will be appreciated that a variety of 
programming languages may be used to implement the teachings of the 
invention as described herein. 

[0060] A machine-readable medium includes any mechanism for storing or 
transmitting information in a form readable by a machine (e.g., a computer). 
For example, a machine-readable medium includes read only memory 
("ROM"); random access memory ("RAM"); magnetic dislc storage media; 
optical storage media; flash memory devices; electrical, optical, acoustical or 
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other form of propagated signals (e.g., carrier waves, infrared signals, digital 
signals, etc.); etc. 

Overview 

[00611 Figure 29 is a block diagram of one emodiment of a coder. Referring 
to Figure 29, a data interface 2901 is coupled to receive data to be encoded or 
output data after decoding. A DC level shifter 2902 is coupled to data 
interface 2901 to perform DC level shifting during encoding and decoding. 
Wavelet transform 2903 is coupled to a DC level shifter to perform forward 
or inverse wavelet transforms depending on the direction of flow. In one 
embodiment, wavelet transform 2903 performs a 5,3 reversible wavelet 
transform and a 5,3 irreversible wavelet transform to decompose an image 
into two to five levels of decomposition. Line buffers 2904 are coupled to 
wavelet transform 2903 to provide memory to support wavelet transform 
2903 when performing the wavelet transform. 

[0062] A scaler quantization/dequantization block 2905 is coupled to 
wavelet transform 2903 to perform scaler quantization. In one embodiment, 
scaler quantization is only used for the 5,3 irreversible wavelet transform. A 
precoder 2906 is coupled to scaler quantizer to perform precoding. In one 
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embodiment, the preceding converts coefficients from two's complement to 
sign magnitude (or vice versa for decoding). The preceding also determines 
zero bitplanes. Work memory A and work memory B are coupled to 
precoder 2906 along with packet header processing 2907. The interface to 
work memory A and B and packet header 2907 is also coupled to bit 
modeling MQ-coders 29081.^. Each of the MQ-coders 2908i.^ is coupled to 
individual code memory 291 1^ to store coded data (compressed data in 
JPEG 2000 terminology). The coded data from code memory and the packet 
header from packet header 2907 is output as coded data. This is a P'EG 2000 
bitstream. An additional functional block (not shown) may be used to 
create/read the main and tile-part headers. The bitstream and the headers 
make a JPEG 2000 codestream. 

Context Model Data Structure with Skipping for Sub-bit-planes 
[0063] In JPEG 2000, for bit-planes of coefficients that are not initially all- 
zero bit-planes (starting with the most significant bit (MSB) down), each 
coefficient is coded in one of three sub-bit-plane passes: significance 
propagation, refinement, and cleanup. Figure 2A shows an example 8x8 
code-block with the sub-bit-plane pass identified for one bit-plane for each 
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coefficient. Referring to Figure 2A, SP refers to the significance propagation 
pass, R refers to refinement pass, and C refers to the cleanup pass. The 
indices from 0 to 63 in Figure 2A show the code block scan order. Thus, the 
scanning order is down four coefficents and then back up to the top row, 
respectively continuing across the code-block. Once scanning has completed 
across the code-block, scanning continues at the fifth coefficient of each 
column down across rest of the code-block. 

[0064] A typical implementation might read the entire block of coefficients 3 
times, once for each coding pass of a bit-plane. This technique describes a 
way to read the entire block of coefficients for the significance propagation 
pass of each bit-plane, but to read only those coefficients actually needed for 
the refinement and cleanup passes. 

[0065] On the left side of each cell in Figure 2A, solid lines indicate 
coefficients in the refinement sub-bit-plane pass and dotted lines indicate 
coefficients that are skipped in the refinement sub-bit-plane pass. Solid lines 
on the right side of each cell in Figure 2A are similar for the cleanup sub-bit- 
plane pass. Once the pass has been identified for each coefficient, then the 
coefficients may be processed. 

[0066] Data structures may be built during the significance propagation pass 
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using a process described below. The data structoes may be used by the 
context model to reduce the number of accesses to the memory. Using the 
data structures, instead of checking each cell to determine what pass the 
irtformation is in and then having to skip them, the context model need only 
access the memory once. Furthermore, the data structures allow access to 
multiple locations at the same time, such as when cleanup bits are coded 4 at 
a time. 

[0067] Tables 1 and 2 show data structures that describe the locations of 
coefficients in the refinement and cleanup sub-bit-plane passes 
(respectively). For each index, there is a run coimt of the number of 
coefficients in the sub-bit-plane pass and a skip count of the number of 
following coefficients that are in a different pass. These data structures 
allow these sub-bit-plane passes to be coded efficiently by allowing 
coefficients in other passes to be skipped. 

Table 1 - Data structure for refinement bits 



index 


run skip 


0 


0 1 


1 


2 8 


2 


1 3 


3 


1 48 
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Table 2 - Data structure for clean.up bits 



index 


run 


skip 


0 


0 


8 


1 


2 


2 


2 


2 


2 


3 


2 


2 


4 


16 


1 


5 


3 


1 


6 


3 


1 


7 


3 


1 


8 


15 


0 



[0068] The process is performed by processing logic that may comprise 
hardware, software, or a combination of both. In one embodiment, the 
processing logic that creates the data structures and uses them is located in 
bit modeling MQ-coders 2908i.^ of Figure 29. During the process, to create 
these data structures, the datastructures are first initialized. 



ri = 0 
ci = 0 

r_run[ri] = 0 
r_skip[ri] = 0 
c_rxm[ci] = 0 
c_skip[ci] = 0 
state ^ INITIAL 
CLEANUP 



/ / index for refinement 

// index for cleanup 

// run coimt for refinement 

// skip count for refinement 

/ / run coimt for cleanup 

/ / skip coimt for cleanup 

// state can be INITIAL, SIG_PROP, REFINE or 



[0069] The "state" variable is used to distinguish between the start and the 
middle of a rim. The "state" variable indicates the coding pass for the 
previous coefficient. If the current coefficient is the same, the size of a run or 
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skip will be increased; if it is different, a new run will be started. 

Then each coefficient in the code-block is considered in the code block scan 

order in order to generate the seperate counts. 

For-yl = 0 to maximum-for-yl step 4 

for X = 0 to maximum-for-x step 1 

for y2 = 0 to maximum-for-y2 step 1 
process coefficient [x,yl+y2] 

In the above, the maximum for y 1 is the greatest integer multiple of 4 that is 

less than the height of the code-block ("(height -1) & -3"). The maximum for 

X is "width-1" for the code-block. The maximum for y2 is the smaller of 3 or 

"height-yl-1". One embodiment of a procedure for processing each 

coefficient is: 

if coefficient was significant in previous bitplane then 
if state is not REFINE then 
ri = ri + 1 
r_nm[ri] = 1 
r_skip[ri] = 0 
state = REFINE 

else 

r_run[ri] = r_nm[ri] + 1 
c_skip[ci] = c_skip[ci] + 1 
else if a neighbor of the coefficient is significant then 
(coefficient is predicted sigificant, code it) 
r_skip[ri] = r„skip[ri] + 1 
c_skip[ci] = c_skip[ci] + 1 
state = SIG_PROP 

else 

if state is not CLEANUP then 
ci = ci + 1 



25 



c_run[ci] = 1 
c_skip[ci] = 0 
state - CLEANUP 

else 

c_nm[ci] = c_rxm[ci] + 1 
r_skip[ri] = r_skip[ri] + 1 

[0070] The result of applying the procedure is that all the coefficients in the 

significance propagation pass will be coded while the data structures for the 

refinement bits and clean up bits are created. 

[00711 If desired, the rxin counts can be prevented from wrapping aroimd 

lines. One embodiment of a process to prevent wrapping aroimd lines is 

described below in the following psuedo-code. This may allow for handling 

boundaries more easily. 

for yl = 0 to maximum for yl step 4 

for X = 0 to maximum for x step 1 

for y2 = 0 to maximum for y2 step 1 
process coefficient [x,yl+y2] 
if state is REFINE then 

ri = ri + 1 

r_run[ri] = 0 

r_skip[ri] = 0 

state = INITIAL 
else if state is CLEANUP then 

ci = ci + 1 

c_rxm[ci] = 0 

c_skip[ci] = 0 

state = INITIAL 

[0072] For software, storing run and skip values as integers (in 32-bits for 32- 
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bit computers, etc.) is most convenient. The worst case is runs of length one, 
with a zero length run at the start. In JPEG 2000, code-blocks are limited to a 
maximum of 4096 coefficents. Also the width and height of code-blocks are 
limited to a maximum of 1024 coefficients. For any size code-block with 

4096 total coefficients with nm counts continuing across groups of lines, 

4097 memory locations is the maximum number of memory locations for the 
memory size. For 64x64 code-blocks with rim coimts starting every group of 
four lines, (4x64+l)x(64/4) = 4112 memory locations is the maximum. For 
1024x4 code-blocks with rim coimts starting every group of four lines 
(4x4+l)x(1024/ 4) = 4352 memory locations is the maximum. 

[0073] In hardware to save memory, a minimal but fixed number of bits can 
be used for run and skip counts. If an indication of whether the first count is 
a run or skip is signaled (e.g., one bit signal indication), then run counts will 
be greater than 1 (and the ability to code 0 is not needed). For any size code- 
block with 4096 total coefficients with run counts continuing across groups 
of lines, one bit may be used to signal whether the first count is a run or skip 
and 4096x12 bits for counts for a total of 49,153 bits. For 64x64 code-blocks 
with run counts starting every group of four lines, one bit may be used to 
signal whether the first count is a run or skip for each group of four lines. 
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Thus, a number of bits would be 1x64/4 + 4096x12 = 49,168 bits. For 1024x4 
code-blocks with run counts starting every group of four lines, the number 
of bits would be 1x1024/4 + 4096x12 = 49,408 bits. 

[0074] One embodiment of a variable length code can be used to represent 
coxmts. Table 3 shows an example where small coimts are represented with 
a small nvimber of bits (e.g., 3 bits) and large counts are represented with a 
larger number of bits (e.g., 13 bits). The goal of such an approach is to have 
most cotmts be either 1, 2, 3, or 4 so that the smaller codewords are used 
more frequently. Only two sizes are used to make the implementation 
simple. However, more than two sizes could be used with added 
complexity. 

Table 3 - An Example of Simple Variable Length Code for Counts 



count 


codeword 


1 


000 


2 


001 


3 


010 


4 


Oil 


5 


1000000000000 


6 


1000000000001 


4096 


1111111111011 



[0075] For this code, the worst case is when all run lengths are 1 (every 
codeword is 3 bits). For three cases (counts go across lines, 64x64 code 
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blocks with groups of four lines, 1024x4 code-blocks with groups of four 
lines), the total number of bits are 12,289 bits; 12,304 bits; and 12,544 bits, 
respectively. 

[0076] Reduced memory usage can be achieved with a more complicated 
variable length code. A good structured code is a gamma code, or / 
(from Appendix A of Bell, Cleary, Whitten "Text Compression", Prenice 
Hall, NJ, 1990.) as shown in Table 4. 



Table 4 - Structured Variable Length Code for Counts 



coimt 


codeword (y' format) 


codeword (y format) 


1 


0 


0 


2 


10_0 


100 


3 


10_1 


110 


4 


110_00 


10100 


5 


110_01 


10110 


6 


110_10 


11100 


7 


110_11 


11110 


8 


1110_000 


1010100 


9 


1110_001 


1010110 


15 


1110_111 


1111110 


16 


11110.0000 


101010100 


32 


111110.00000 


10101010100 


64 


111110_000000 


1010101010100 


128 


11111110_0000000 


101010101010100 


256 


111111110_00000000 


10101010101010100 


512 


1111111110_000000000 


1010101010101010100 


1024 


11111111110.000000000 


101010101010101010100 


2048 


lllllllllllO.OOOOOOOOOO 


10101010101010101010100 


4096 


llllllllllllO.OOOOOOOOOOO 


1010101010101010101010100 
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[0077] Note that y' and Y only differ in how the bits in the codeword are 
arranged. The "_" in codewords is not part of the codeword, it is just to 
make them easier to read by separating the prefix from the counter. The 
worst case is for coxmts of 2 which requires 3 bits. For the three cases 
(counts go across lines, 64x64 code-blocks with groups of 4 lines, and 1024x4 
code-blocks with groups of four lines), the total number of bits are 6,145 bits; 
6,160 bits; and 6,400 bits, are required, respectively. 
[0078] Table 5 is one embodiment of a code for coimts 1 . . .4096 where the 
longest codeword is 16-bits. Coimts of 1, 2, 3, 4 or 5 are represented with 1, 
2, 4, 5 and 5 bits respectively. Coimts of 6... 21 are represented with 8 bits. 
Counts of 22. . .4096 are represented with 16-bits. Coimts of 3 and 6 are the 
worst case. For the three cases (counts go across lines, 64x64 code-blocks 
with groups of 4 lines, and 1024x4 code-blocks with groups of 4 lines), the 
total number of bits are 5,463 bits; 5,478 bits; and 5,718 bits. 
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Table 5 - Optimal Variable Length Code for Coiints 



count 


codeword 


1 
1 


u 


2 


01 


3 


1100 


4 


1101_0 


5 


1101_1 


6 


1110_0000 


7 


1110_0001 


21 


1110_1111 


22 


1111_0000_0000_0000 


23 


1111_0000_0000_0001 


4096 


1111_1111_1110_1010 



[0079] If using variable length codes with hardware, accessing both a run 
count and a skip count during the same clock cycle may be desirable. In one 
embodiment, the memory is larger than the minimum size by one word 
(e.g., 16 bits) so that flushing is simple. For example, with 16 bit words, run 
coimts use 163 bits. Therefore, the last word uses only 3 bits. This requires 
padding with 13 bits to complete a word. Skip counts use 85 bits. Therfore, 
the last word uses only 5 bits. This requires padding with 11 bits to 
complete a word. If 



memory size > 



163 + 13 85 + 11 
- + - 



.16 16 ; 

or 17, padding of rim and skip coimts can be independent. If the size is 16, 
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there is one word with both run counts and skip counts, with padding in 
between. Figure 30 illustrates an example of a 16 bit word having both rim 
counts and skip counts. 

[0080] Figure 2B illustrates a memory for a variable length run and skip 
coimts. This memory structure allows starting a run coimt from one side 

(201) of memory and skip counts from the other side (202) of memory. One 
side 201 of memory may be the beginning or the end, with the other side 

(202) being the end or beginning, respectively. This allows for 
simultaneously decoding a rim count and a skip count in parallel because 
the start of both 1 run and 1 skip variable word is known at the same time, 
instead of decoding a nm coimt first, determining its length, and then 
decoding a skip covint. If they were serial (run-skip-run, etc.), one at a time 
would be known (as one would have to be decoded to find the other). 
[0081] Note that separate skip coxmt and nm count decoders may be used to 
decode the information in the memory shown in Figure 2B. Alternatively, a 
single decoder could be shared by both. 

Hardware for Context Model 

[0082] The context model may be implemented in hardware. In hardware, 
one goal is to generate the next context as soon as possible so the MQ coder 
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is not idle. 

Memory Organization 

[0083] Figures 3A-D illustrate the neighborhood coefficients and memory 
organization for embodiments of a context model. The context model for 
individual coefficients is at most based on a 3x3 neighborhood as shown in 
Figure 3A. In one embodiment, four bits are handled at the same time. In 
such a case, the context model for groups of four coefficients 311-314 is at 
most based on a 3x6 neighborhood as shown in Figure 3B. It is often 
preferable for memory accesses to hardware to be dealing with data 
grouped in powers of two. Therefore, the region that is based on a power of 
two that would contain the 3x6 region is a 4x8 region. Figure 3C illustrates a 
3x6 neighborhood being a subset of a 4x8 region of coefficients. Access to 
the entire 4x8 region of Figure 3C may be performed as separate accesses 
that may or may not be to the same memory. Figure 3D illustrates a 4x8 
region that is divided into four 2x4 regions 301-304 which are each stored in 
a different memories for parallel, random access. This memory structure 
allows everything needed to determine contents from a group of coefficients 
to be read from memory at once, instead of sequentially. That is, the entire 
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4x8 block of coefficient information may be accessed at once. 
[0084] Figure 4 shows one embodiment of the significance memory 
organization for random access for a 16x16 code-block. One implementation 
may handle other sizes, including, but not limited to, 32x32 and 64x64. 
Referring to Figure 4, each coefficient has been assigned to one of four 
memories (A, B, C or D). Certain groups (two rows along the top and two 
rows along the bottom) are half the size of the other groups. This is because 
in Figure 3D, the top two rows are outside of the code-block (off the edge) 
for the first row of code-blocks. There is a similar boimdary condition of the 
bottom of the code-blocks. 

[0085] In one embodiment, these memories store one bit per coefficient for 
significance state (8 bits total per address). In an alternative embodiment, 
these memories store two bits per coefficient for significance state and sign 
(16 bits total per address). In still another embodiment, these memories 
store the entire coefficient (8N bits if N is the size of one coefficient). In an 
alternative embodiment, if the entire coefficient is not stored in these 
memories, an additional single memory with one address per coefficient 
may be used. 

[0086] The following Verilog code is one embodiment of a procedure to 
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convert two 6-bit address inputs ("x" and "y") from control logic in the 

context model (for the 64x64 code-blocks) into a 7-bit address output to 

memory ("addr") and a 2-bit "bank" select to indicate memory "A", "B", 

"C" or "D". 

module makeAddress(x, y, addr, barJ<:); 

input [5:0] x; /*x has bits 5. . .0, where 5 is MSB 

input [5:0] y; 
output [6:0] addr; 
output [1:0] bank; 

wire [5:0] yp2; 

assign yp2 = y + 2; 
assign addr = {yp2[5:3], x[5:2]}; 
assign barvk = {yp2[2], x[l]}; 
endmodle 

[0087] The first assign statement above sets the offset for the boundary. In 
other words, the offset "assign yp2 = y + 2" is used for proper aligimient of 
groups of four as shown in Figure 4. The second assign statement above sets 
the address to bits 5..3 of the input y plus the offset concatenated with bits 
5..2 of input x as the lower poriton of the nimiber. The third assign 
statement above sets the bank equal to the bit 2 of the input y plus the offset 
concatenated with bit 1 of input x. 
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Significance Propagation Pass 

[0088] Figure 5 shows one embodiment of the memories and registers used 
in the significance propagation pass for random access. Referring to Figure 
5, address A is input to memory A to produce data output, which is also 
stored in register 501. In response to address B, memory B outputs data 
which is also stored in register 502. Similarly, memory C outputs data in 
response to address C, where the output is also stored in register 503. 
Lastly, memory D outputs data in response to address D and also stores the 
data in register 504. In one embodiment, the outputs of each of the 
memories A-D is a 2x4 region which together create a 4x8 region (e.g., region 
601 of Figure 6). 

[0089] All of the outputs of the memories and registers in Figure 5 together 
provide a 6x6 region of significance bits. Note this could be significance 
state and sign or the actual coefficients in alternative embodiments. In other 
words, the data being used out of the memories A-D in parallel is combined 
with the data read outs from the memories A-D in the previous cycle that 
had been stored in registers 501-504. This region of significance bits plus 
feedbaclc from the context model is enough to determine in which pass a 4x4 
region of coefficients is. 
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[0090] Figure 6 illustrates the significance state from memory and stored in 
registers for random access. Referring to Figure 6, region 601 illustrates a 
4x8 region that is read from memory. Region 602 illustrates a 3x6 region 
that is read from memories A-D and is used for context modeling. Region 
603 illustrates a 3x6 region that is stored in registers 501-504 and is used in 
context modeling. Region 604 illustrates a 4x4 region of coefficients that are 
being processed. Figure 6 also shows exemplary 2x4 portions of the 8x8 
block that may come from memory locations of memories A-D and registers 
501-504. 

[0091] One embodiment of address generation logic for the significance 
propagation pass is set forth in the following pseudo code. Note that the 
addressing is not dependent on the data, and zero data is provided at 
boimdaries. 

address_A_y = 0 

address_B_y = 0 
address_C_y = 0 
address_D_y = 0 
f or y = 0 to 60 step 4 

address_A_x = 0 

address_C_x = 0 

read memory A (will be registered next) 

read memory C (will be registered next) 

assert clear for memory B register (will be cleared next) 

assert clear for memory D register (will be cleared next) 
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for X = 0 to 60 step 4 

address_A_x = x+4 

address_B_x = x 

address„C_x = x+4 

address__D_x = x 

if X < 60 then 

read memory A (will be registered next) 
read memory C (will be registered next) 

else 

use "all bits zero" for memory A output 
use "all bits zero" for memory B output 
read memory B(will be registered next) 
read memory D(will be registered next) 
process 4x4 block of coeffcients x. . .x+3, y . . .y+3 
if y AND 4 ==0 

address_A_y = address__A_y + 8 
address__B_y = address_B_y + 8 

else 

address_C_y = address_C„y + 8 
address_D_y = address_D_y + 8 

[0092] For processing 4x4 blocks, runs of bits in the same pass may be 

handled together. If there are N coefficients in a row in the refinement pass, 

the following pseudo code can be used to process them: 

if state is not REFINE then 
ri = ri + 1 
r_run[ri] = N 
r_skip[ri] = 0 
state = REFINE 

else 

r_nm[ri] = r„rim[ri] + N 



Note that this code is similar to the code described above except "N" has 
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replaced 1 to signify that N coefficients are being processed instead of 1. 
[0093] If there are N coefficients in a row in the cleanup pass, the following 
pseudo code depicts one embodiment of a process for processing 
coefficients: 

if state is not CLEANUP then 
ci = ci + 1 
c„run[ci] = N 
c_skip[ci] = 0 
state - CLEANUP 

else 

c_run[ci] - c_run[ci] + N 
[0094] Figure 7 is a block diagram of one embodiment of significance 
propagation pass logic. In one embodiment, this logic is in bit modeling 
MQ-coders 2908i.^ of Figure 29. The pass for each coefficient is for a 4x4 
region represented as three bits per coefficient: significance propagation or 
other, refinement or other, and cleanup or other. By controlling the access 
operation of the memory A-D in Figure 5, a 4x4 block is retrieved from 
memory and the significance propagation pass is performed. Looking at 
each 4x4 block, the runs in the various passes are identified and any 
coefficient in the significance pass is to be coded, while rims and skip counts 
for the refinement and clean up passes are identified for processing one rim 
at a time. When the previous bit in the block (via the scan order) (or a 
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previous codeblock if starting a new codeblock) is a significance propagation 
bit and the current state is not significance propagation, a new run is started. 
In such a case, the index is incremented in the table storing run and skip 
counts (e.g., set skip to zero, set run to the first value). Both tables of run 
and skip counts are incremented in this manner while processing the 4x4 
block. If the previous bit in the 4x4 block was in the refinement or cleanup 
passes and more such data foUows, then the count of the current rvin is 
incremented. Other sized regions may be used, including, for example, a 
4xN region, where N is a power of two. 

[0095] Referring to Figure 7, a significant state 701 for an 8x8 region is input 
to logic 702 which determines the pass. Figure 31 illustrates an exemplary 
8x8 region. The significance state 701 may comprise information that 
indicates, for example, there are N coefficients in a row in a refinement pass. 
Such information may be accessed from a table, as described above. 
Determine pass logic 702 looks at the 16 3x3 regions in the center 6x6 region 
of the 8x8 region. Coefficients A..I represent the first 3x3 region. Figure 32 
is a diagram of one embodiment of determine pass logic 702. The logic in 
Figure 32 is replicated 16 times, one for each coefficient in the 4x4 block. 
Note that the regioris may be other than 3x3 regions and the number of 
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regions processed may be greater than or less than 16 at a time. 
[0096] Referring to Figure 32, all bits of coefficients A-C are input to OR gate 
3201. All bits of coefficients D and F are coupled to inputs of OR gate 3202. 
All bits of coefficients G-I are coupled to inputs of OR gate 3203. The 
outputs of OR gates 3201-3203 are coupled to inputs of OR gate 3204. The 
output of OR gate 3204 is coupled to an input of inverter 3206 and an input 
of AND gate 3208. Coefficient E represents the 16 bit output of refinement 
signals 704 and is coupled to an input of inverter 3205, the output of which is 
coupled to the other input of AND gate 3208 and an input of AND gate 3207. 
The output of AND gate 3208 is the significance propagation signals 703. 
The output of inverter 3206 is coupled to the other input of AND gate 3207. 
The output of AND gate 3207 is the cleanup signals 705. 
[0097] In operation, if any of the significance state bits E are 0, then the 
output of AND gates 3208 corresponding to those bit positions, and thus, the 
significance propagation signals 704, are going to be I's if the significance 
state is 1 for any of coefficients A-D or F-L Similarly, if any of the 
significance state bits E are 0, then the output of AND gate 3207 
corresponding those bit positions, and thus the cleanup signals 705, are 
going to be I's if the significance state bits are all zero for coefficients A-D or 



41 

F-I. 

[0098] As a result of the determination, logic 702 asserts either significant 
propagation signals 703, refinement pass signals 704, or clean up pass 
signals 705. In one embodiment, each of signals 703-705 is 16 bits wide. For 
each corresponding bit in signals 703, 704, and 705, only one bit may be 1 
and 2 must be 0. Thus, there are three possibilities at 16 locations. Each of 
the outputs of logic 702 is coupled to one input of selection logic (e.g., 
multiplexer (MUX)) 707. 

[0099] Selection logic 707 generates three pass bits for the current coefficient 
indicating the pass for the current coefficient and sends the pass bits to 
control logic 709. Only one of these pass bits is asserted in response to count 
signal 708 output from control 709. Coimt signal 708 indicates which of the 
16 coefficients in the 4x4 block is currently by processed. When handling 
refinement bit rtms and cleanup bit rtms, coimt signal 708 may be 
incremented by a number greater than one. Thus, of the 16 bits in each of 
the outputs of determine pass logic 702, the bit in each of the three outputs 
corresponding to that coefficient is output. 

[0100] The refinement pass signals 704 and clean up pass signals 705 are 
input to mask 705 with feedback coimt signal 708. Count signal 708 is the 
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current coefficient position in the 4x4 region, e.g., 0. ..15. In response to 
these inputs, mask 705 masks what has already been done, as indicated by 
count 708, including only coefficients that are not being coded yet. For 
example, if three coefficients have already been processed, mask 705 masks 
the three signal lines of each of the refinement and cleanup outputs (704 and 
705). 

[0101] Mask 705 also generates 2 outputs to priority encoder 706 
representing signal 704 and 705 with certain signals masked to one (in one 
embodiment). These outputs of mask 705 are masked refinement and 
masked cleanup indications (e.g., signals). 

[0102] In response to the two inputs, priority encoder 706 finds the next non- 
refinement bit (or coefficient) and the next non-cleanup bit for the 
significance propagation pass and inputs these to control logic 709. In one 
embodiment, priority encoder 706 is a zero-finding priority encoder. In so 
doing, priority encoder 706 converts the current position of the bit (or 
coefficient) in the codeblock into a count of leading zeros. In one 
embodiment, this is performed using a truth table such as 
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Input 


Output 


1 


X X X X 


X 


0 


0 


1 X X X 


X 


1 


0 


0 1 X X 


X 


2 



[0103] Mask 705, priority encoder 706, and selection logic 707 may comprise 

a processing tmit that receives outputs from determine pass unit 702 and 

generates outputs indicating the next non-refinement coefficient and next 

non-cleanup coefficient and the pass for the current coefficient. 

[0104] In response to its inputs, control logic 709 generates a refinement next 

index, a refinement rim indication, a refinement slcip indication, a cleanup 

next index, a cleanup rim indication, a cleanup skip indication, and a 

significant propagation indication. The inputs to control logic are: 

next non-refinement bit position "R'' 
next non-cleanup bit position ''C" 

If R > coimt then 

refinement run = R-coimt 
refinement skip = 0 
clean-up rim = 0 
clean-up skip = R-cotint 
refinement next index =1 
cleanup next index =0 
signif prop =0 

Else if C > count 

refinement run = 0 
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refinement skip = C-count 
clean-up run = C-count 
clean-up skip = 0 
refinement next index = 0 
clean-up next index =1 
signif prop = 0 

refinement rim = 0 
refinement skip = 1 
clean-up run = 0 
clean-up skip = 1 
refinement next index = 0 
clean-up next index =0 
signif prop = 1 

[0105] The following pseudo code describes the operation of the significance 

propagation pass logic described in Figure 7. 

coimt = 0 

while (coimt < 16) 

mask = (1 « count)-l 

refinement_masked = refinement I mask 

use priority encoder to find next non-refinement bit 

cleanup_mask = clean_up I mask 

use priority encoder to find next non-cleanup bit 

if current bit is in significance propagation pass 

process coefficient as significance propagation 

coimt = count + 1 
else if current bit in refinement pass 

N = "next non-refinement bit'' - count 

process N bits as refinement pass 

count = coxmt + N 

else 

N = "next non-cleanup bit" - count 



else 
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process N bits as cleanup pass 
count = count + N 

[0106] Note that the significance states are updated from the MQ decoder 
(and MQ encoder or coefficient values when encoding) anytime a "1" 
coefficient is coded in the significance propagation pass. 
[0107] Assuming the context model operates in one clock cycle and the MQ 
coder operates in one clod?; cycle, two clock cycles are required when there is 
feedback. Figure 8 illustrates an example of performance on a 4x4 block, 
illustrating the potentially worst case scenario. Eight context models and 
MQ coders operating in parallel at twice the component clock rate should be 
able to decode 7 bitplanes per coefficient (8x2/2.25 « 7). Without skipping in 
the significance propagation pass, the worst case performance reduces to at 
most 5.5 bitplanes per coefficient. Without skipping in any pass, the worst 
case performance reduces to at most 4 bitplanes per coefficient. 

Significance Propagation Pass Skipping in Software 

[0108] In software, parallel access from many memories is impractical. 
Instead, in one embodiment, the code-block is divided into 4x4 groups of 
coefficients. For each group, a coimt is maintained of the number of bits that 
are significant. In such a case, the maximum memory required is 256x5 bits. 
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Blocks of coefficients that are all in the refinement pass have a count of 16. 
Blocks that have coimt 0 may be all cleanup and just need to have their 
neighbors checked to see if they are all cleanup. 



Cleanup and Refinement Passes 

[0109] For the cleanup pass, addressing is data dependent and may be 

generated using the following psuedo code. The address x,y of the next 

coefficient in the cleanup pass is input. 

module cleanup Address(x, y, addrA, addrB, addrC, addrD) 
input [5:0] x; 
input [5:0] y; 
output [6:0] addrA; 
output [6:0] addrB; 
output [6:0] addrC; 
output [6:0] addrD; 

wire [5:0] yp2; 
wire [4:0] ax; 
wire [4:0] bx; 
wire [4:0] cx; 
wire [4:0] dx; 

assign yp2 = y+2; 

assign ax = (x[l:0] == 3) ? x[5:2l + 1 : x[5:2]; 
assign cx = (x[l:0] == 3) ? x[5:2] + 1 : x[5:2]; 
assign bx = (x[l:0] == 0) ? x[5:2] - 1 : x[5:2]; 
assign dx - (x[l:0] == 0) ? x[5:2] - 1 : x[5:2]; 

assign ay = y[2] ? yp2[5:3] + 1 : yp2[5:3]; 
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assign by = y[2] ? yp2[5:3] + 1 : yp2[5:3]; 
assign cy = yp2[5:3]; 
assign dy = yp2[5:3]; 

assign addrA = {ay, ax}; 
assign addrB = {by, bx}; 
assign addrC = {cy, cx}; 
assign addrD = {dy, dx}; 
endmodule 

[0110] The addressing used for the cleanup pass may also be used for the 

refinement pass. However, a smaller neighborhood is sufficient for the 

refinement pass. 

If (yp2[l:0] == 1) or (yp2[l:0] == 2) then 
ifyp2[2] ==lthen 

just read memories C and D 

else 

just read memories A and B 

else 

read memories A, B, C and D 



Sequential Addressing for All Passes 

[0111] With sequential addressing for all passes, a simpler memory 
organization using two memories can be used. Figure 9 illustrates one 
embodiment of an organization of a significance memory for sequential 
accesses of 4x4 regions of a code-block. Referring to Figure 9, each 4x4 
region is assigned to one of two memories A or B. This allows parallel 
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access for everything needed for a 16x16 block. Note that the first code- 
blocks are only half because the offset is similar to that described above and 
only two rows of coefficients are relevant when processing 8x8 blocks such 
as in Figure 6 where the top two rows do not include actual data being 
processed. 

[01121 Figure 10 illustrates one embodiment of the memories and registers of 
a memory path used for the significance propagation pass. Referring to 
Figure 10, memory A in response to address A generates data output. 
Similarly, memory B generates data output in response to address B. A 2x2 
crossbar 1003 has inputs coupled to the outputs of memories A and B. One 
output of crossbar is coupled to register 1001 and one output of the memory 
path. The other output of crossbar 1003 is coupled to register 1002 and the 
other output of the memory path. Thus, the outputs of memories A and B 
may be stored in either registers 1001 and 1002 and on either output. Data 
read out from memories A and B is for a 4x4 region. Registers 1001 and 1002 
store a 5x4 region. When registers are loaded, the rightmost 1x4 column is 
moved to the leftmost 1x4 column and the other columns are loaded from 
the memory data output. Crossbar 1003 controls the output of the data from 
memories A and B to the proper outputs of the memory path via "ping- 
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ponging'' data into the outputs as data is processed row after row. 
[0113] Figure 11 shows how memory and registers of Figure 10 are used to 
provide the proper region for context model operation. Referring to Figure 
11, region 1102 is a 4x4 region of coefficients to be processed. Region 1101 
represents a 5x6 region stored in registers 1001 and 1002 used for context 
modeling (with the 5x1 regions above and below the 5x6 region being 
ignored). Region 1103 is a 4x8 region from memory. Region 1104 is a 1x6 
region from memory used for context modeling. 

[0114] One embodiment of the pseudo code for addressing memory for all 
three coding passes is as follows. 

address_A_y = 0 
address_B_y = 0 
for y = 0 to 60 step 4 

address_A_x = 0 

address_B_x = 0 

clear registers 

read memory A (will be registered next) 
read memory B (will be registered next) 

for X = 0 to 60 step 4 

address_A_x = x+4 

address_B_x = x+4 

if X < 60 then 

read memory A (will be registered next) 
read memory B (will be registered next) 

else 

use "all bits zero" for memory A output 
use "all bits zero" for memory B output 
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process 4x4 block of coeffcients x. . .x+3, y . . .y+3 
if y AND 4 ==0 

address_A_y = address_A_y + 8 

else 

address_B_y = address_B_y + 8 
[01151 The memory contams state to mdicate the correct pass for the 
refinement and cleanup passes. The state may be 2 bits per coefficient to 
distinguish between the three states (significance propagation, clean-up and 
refinement). 

[01161 During the significance propagation pass, in parallel for all 16- 
coefficients, the state is set to refinement for all significant coefficients and to 
cleanup for all non-significant coefficients. As the processing on the 16- 
coefficients continues, the state of any coefficients that are in the significance 
propagation pass is changed from cleanup to significance propagation. The 
state may be 1 bit per coefficient, referred to herein as the "pass" bit. In one 
embodiment, the significance state and the pass bit are used to determine the 
correct pass. Table 6 below illustrates use of the pass bit. Since 1 bit per 
coefficient is used, this uses less memory than the rtin and skip count 
methods described herein. 
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Table 6 - Using "Pass" Bit 



pass 


significance state, comment 
pass bit 
current next 


significance 
propagation 


l,x 1,1 in refinement pass* 

in parallel for all 16 coefficients at 
start 

0,0 0,0 in cleanup pass 

0,1 code "0" in significance propagation 
1,0 code "1" in significance propagation 

0,1 0,1 code "0" in significance propagation 
1,0 code "1" in significance propagation 


refinement pass 


14 14 


cleanup pass 


0,0 0,0 still not significant 
l,x became significant 



*In one embodiment, this is done in a parallel for all 16 coefficients at the 
start of processing for the 4x4 block. 

[0117] The memories providing access to 2x4 regions can be 48 bits wide 
with sigruficance state, pass and sign bits for each coefficient. 
[0118] Figure 12 is a block diagram of one embodiment of pass 
determination logic that implements Table 6 above and uses a priority 
encoder to find each coefficient that is in the current pass. Referring to 
Figure 12, determine pass logic 1203 receives a significant state 1201 for a 6x6 
region, a pass bit 1202 for a 4x4 region, and a current pass signal (or other 
indication) 1220 that indicates the current pass. Pass bit 1202 comprises a 
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signal for each coefficient in the 4x4 region (i.e., 16 signals). In response to 
these inputs, determine pass logic 1203 generates an output to indicate the 
pass for the 4x4 region. In so doing, for each coefficient in the 4x4 region, 
determine pass logic 1203 asserts signals 1204 to indicate a sigruficance 
propagation pass bit for coefficients in the significance propagation pass, 
signals 1205 to indicate a refinement pass bit for coefficients in the 
refinement pass, or signals 1206 to indicate a cleanup pass bit for coefficients 
in the cleanup pass. 

[0119] Selective logic 1207 is responsive to current pass signal 1220 to output 
one of the indications 1204-1206 to mask logic 1208. In one embodiment, 
selective logic 1207 comprises a 16x3:1 multiplexer (Mux). Mask logic 1208 
generates a signal in response to coimt signal 1210, which indicates the 
coefficient currently being processed. The output of the mask 1208 is input 
to priority encoder 1209 which outputs the signal to control logic 1212. 
Mask logic 1208 and priority encoder 1209 operates in a maimer similar to 
their named coxmterparts in Figure 7. In response to the inputs, control logic 
1212 generates a code or idle indication on signal line 1213 and a coimt 
signal 1210. 

[0120] Next pass bit logic 1211 receives an output from priority encoder 1209 
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(indicating the position of the coefficient currently being processed), the 
current pass signals 1220, a new significant state 1221 from an MQ coder, 
and the refinement pass signals 1205. The refinement pass signals 1205 
represent significance state information by indicating whether coefficients 
are previously signficant. Together, current pass signal 1220 and new 
significance state 1221 indicate whether processing is in the cleanup pass. In 
response to the inputs, next pass logic 1211 generates a next pass bit, which 
is used as an output to distinguish between 'code ''0'" and 'code "1"' cases in 
Table 6. The next pass bit may be stored in memory and used subsequently 
as the pass bit 1202. 

[0121] The operation of the logic in Figure 12 is described in the pseudo code 
below. Such functionality may be included in MQ-coders2908i_j^. The 
significance state and pass bit are cleared before processing the first cleanup 
pass. 

count = 0 

if significance propagation pass then 

set next pass bit to "1" for all coeffients in refinement pass 

while (count < 16) 

if significance propagtion pass then 

in_pass = coefficients in significance propagtion pass 
else if refinement pass then 

in_pass = coefficients in refinement pass 
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else 

in_pass = coefficients in cleanup pass 
mask = (1 « count)-l 

in_pass_masked = in_pass AND (NOT mask) 
use priority encoder to find next coefficient in pass, N 
if next coefficient found 
code coeff N 

if significance propagation pass then 

next pass bit = NOT next significance state 
else next pass bit = pass bit 
count = N+1 

else 

count = 16 
idle for coeff N 
next pass bit = pass bit 

[0122] In the above code, "in-pass" variable is the output of a 3:1 

multiplexing function. The "mask" variable represents the mask, and the 

"in-pass-masked" variable represents the resuh of applying the mask. The 

variable N represents the next coefficient in the pass and is the output of the 

priority encoder. Once this has been foimd, the control function in the code 

follows. 

[0123] In the above code "code coeff N" refers to coding coefficient when the 
coefficient is in the current pass. The "idle for coeff N" is executed when 
processing a rim. 



Double Context Generation 



55 

[0124] Contexts often depend on the last bit coded, causing a time critical 
feedback loop between the MQ-coder and context model when decoding. 
To reduce this context model delay to the time for a simple multiplex, the 
context model can generate two contexts, one for when the current bit being 
decoded is "0" and another when it is "1". When the bit is known, selection 
of the context may occur. 

[0125] Figure 13 is a block diagram of one embodiment of double context 
generation logic. Referring to Figure 13, context model generates a context 
0, enable 0, context 1, and enable 1. The context model generates 2 contexts, 
0 and 1, which are both sent to multiplexer (mux) 1302. Mux 1302 is coupled 
to receive these signals and generates a context indication and an enable 
indication to indicate if the context is valid and if a bit should be coded. 
These outputs are coupled to inputs of MQ coder 1303 which generates a bit. 
The output bit from MQ coder 1303 is used by mux 1302 to select the context 
that is to be output to MQ coder 1303. Thus, the context model generates 2 
contexts, one for when the current bit is decoded as "0" and another for 
when it is decoded as "1", and the output bit for the MQ coder 1303 selects 
the correct one. 

[0126] For rim-length coding for the cleanup passes. Table 7 shows the two 
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possible next contexts for each case and follows the cleanup pass flow 
described in the JPEG 2000 Standard. The value of the bit coded in the run- 
length context is used to decide if the next context is for the next group of 
four coefficients or is in the uniform context for the current group of four 
coefficients. If this bit is "0", then the run-length coding succeeded in 
representing four coefficients with one decision and having the MQ-coder 
idle for the following cycle does not significantly hurt speed performance. 
After the bit coded in the second uniform context ("uniform B") is a sign bit 
that is always coded directly in context 9 (with XOR bit 0) in the JPEG 2000 
Standard, where the XOR bit being 0 indicates that the sign is not being 
flipped. 

Table 7 - Double Context Generation for Rim-length Coding in Cleanup 

Passes 



current context 
(being used by MQ 
coder) 


next context for "0" 


next context for "1'' 


run-length 


run-length or 


uniform A 




magnitude 






(or idle) 




uniform A 


imiform B 


xmiform B 


uniform B 


sign 


sign 



A and B above refer to the two bits that are described in Section D.3.4 as 
"[t]he next two bits, returned with the UNIFORM context" in the JPEG 2000 
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standard. 

[0127] For the sigr\ificance propagation and cleanup coding passes, the cases 
without run-length coding are show in Table 8. While a magnitude bit is 
being coded, the contexts are generated for the magnitude of the next 
coefficient assuming the current coefficient is "0" or the sign bit context for 
the current coefficient. 



Table 8 - Double context generation for significance propagation and 

cleanup coding passes. 



current context 
(being used by MQ 
coder) 


next context for "0'' 


next context for "1" 


magnitude bit 


magrutude next 


sign for current 


coefficient 


coefficient 


sign bit 


magnitude next 


magnitude next 


coefficient 


coefficient 



[0128] For the refinement pass, the value of any previous refinement 
coefficient coded does not effect the context. 



MO-Coder 

MQ-Decoder Dataflow with Late Context 

[0129] Figure 14B is a bloclc diagram of a typcial decoding implementation. 
Referring to Figure 14B, context model 1430 provides a context to a memory 
1431 where a probability state is determined. The probability state is 
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converted with logic 1432 to a "Qe_value" for the arithmetic coder 1433, 
which updates an internal A «& C register and determines a decision (MPS or 
LPS). All of this must typically happen before the next context can be 
determined. In many hardware implementations, the decode speed is 
limited by a large feedback loop (feeding back to context model 1431). 
[0130] In contrast, Figure 14A is a block diagram of one embodiment of an 
"early context" MQ-decoder. In this case, the feedback loop has much 
simpler logic 1407, instead of an entire decode operation. Therefore, much 
of the decode and update can be done in parallel with the lower feedback 
loop 1401. 

[0131] Referring to Figure 14A, codestream 1400 is input and updates the 
internal state 1401. In one embodiment, the A and C register of the internal 
state specify a current interval as set forth in the JPEG 2000 Standard in 
Appendix C. The register A indicates the current interval and the code 
register C is the concatination of the Chigh and Clow registers. 
[0132] Context 1402 is provided by context model 1410. Context 1402 is 
used to look up probability state 1404 in memory 1403, which is then 
converted by logic 1405 to a probability class (Qe_value) 1406. Qe_value 
1406 represents the current estimate of a less probable symbol (LPS). 
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Qe_value 1406 is compared with A and C register values as set forth in 
Figure C-15 of the JPEG 2000 Standard of the MQ-coder's internal state by 
logic 1407 to generate the output decision 1408, which may be a more 
probable symbol (MPS) or LPS. The output decision 1408 is input to context 
model 1410. In one embodiment, the operations on the Qe_value and the 
internal state require 16-bit arithmetic. Operations of these blocks 
implement the decoding of a decision as set forth in section C.3.2 of the JPEG 
2000 Standard. 

[0133] Figure 15 is a block diagram of one embodiment of a "late context" 
MQ-decoder. Referring to Figure 15, 16-bit processing has been eliminated 
from the context model feedback loop. Codestream 1501 is received as an 
input to update logic 1503, which updates the internal state 1503, including 
the A and C registers that specify the current interval. The new A and C 
register values and the codestream are input to logic 1504 which generates 
two pclasses as described below, pclass 1509 and pclass 1510, which are 
input to comparison logic 1511 and 1512. 

[0134] Context model 1520 generates context 1502. Context 1502 is used to 
look up a probability state 1506 of memory 1505. In one embodiment, 
memory 1505 comprises a look up table. Identification of probability state 
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1506 allows determining the Qe_Value. Probability state 1506 output from 
memory 1505 is converted by logic 1507 into probability class (index) 1508. 
[0135] Comparion logic 1511 determines if pclass 1509 is greater than the 
probability class index 1508 and comparison logic 1512 compares to 
determine whether probability class index 1508 is greater than pclass 1510. 
The results of both comparison logics 1511 and 1512 are input to AND gate 
1513 such that if both comparisons are true, a decision is output. This 
decision may be an MPS or LPS. Thus, context 1502 is converted into a 5 bit 
probability class index 1508 (since there are 32 possible values for Qe_Value 
in JPEG 2000). The internal state is used to generate two 5-bit probability 
class indices. If the index corresponding to the context is outside the two 
indices generated from the state, the decision is a MPS; otherwise, it is a LPS 
(i.e., inside the two indices). 

[0136] An important advantage of the embodiment of Figure 15 is that the 
internal state update is parallel with generating the next probability class 
(indices) 1508, instead of serially as shown in Figure 14B. Also, because two 
probability classes are only 5 bits when compared to a pclass index, the 
arithmatic is much simpler. 

[0137] Logic 1504 of Figure 15 creates the iiiformation that is depicted in 
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Figure 16A. Given the values in the A and C registers, logic 1504 determines 
what the two split points are for the pclass and then determines if the code is 
in between or outside the split points. These may be done in parallel. 
[01381 Figure 16A illustrates how the comparison of probability class indices 
works. Referring to Figure 16A, "pclass 0" is a high skew case with most of 
the interval devoted to MPS. For "pclass 1" though "pclass 4", the skew is 
less and the MPS interval shrinks, "pclass 5" shows the MPS a conditional 
exchange which occurs for probabilities close to 50%. The "known state" has 
a codestream value ("code") which would be a MPS for some probability 
classes and a LPS for others. Because the probability classes are ordered, 
two comparisons are sufficient to determine if "code" is a MPS or LPS. In 
other words, in Figure 16A, given the location of the code in the known 
state, the decision will be an MPS for pclass 0-3, but will be an LPS always 
for pclass 4, and then again an MPS for pclass 5. Instead of figuring out 
whether it will be an MPS or LPS for each probability class, only the two 
breakpoints (between pclass 3 and 4 and between pclass 4 and 5) need to be 
determined. Therefore, when the QE value is given (when ihe probability 
class/ index is known), a determination may be made as to what probability 
class is actually in the space where the break points are. 
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[0139] A similar method in hardware could be used to determine the MPS or 
LPS for each possible Qe_Value and then multiplex the result. For example. 
Figure 16B shows a mux 1610 having a number of inputs, each of which are 
associated with a pclass and provides as an output either an MPS or LPS 
depending on the code. 

Multiple Bit Decoding with MQ-coder 

[0140] Multiple MPSs can be decoded (for consecutive uses of the same 
PClass) at one time as long as none of the MPSs or only the last MPS 
requires normalization. Figure 17 illustrates the intervals for multiple MPS 
decoding. In the standard, if the difference between where the codestream is 
in relation to the interval specified by the A and C registers and the 
Qe_value is 2 or greater, then multiple MPSs can be decoded. If the interval 
size is divided by Qe_value and then if the decoder remains in the same 
context and, thus, the same probability class, then multiple MPSs can be 
decoded at one time. For example, when looking at the codestream and 
knowing the 16 bits that are being processed, if the location of the codesteam 
in the interval specified by the A and C registers is multiple Qe_values 
away, indicating that the same context is going to be used to process the 
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data for multiple cycles, and thus the same probability class, then multiple 
MPSs may be decoded at one time. In other words, if 

inter val specified by A and C registers - location of the codestream 

is determined and roimded to the next lowest interger, the result indicates 
the number of MPSs that may be decoded at one time. This calculation may 
be performed by well-known hardware. 

An Exemplary Implementation of a 5.3 Filter 

[0141] In one embodiment, a reversible and irreversible 5,3 wavelet filters 
are used. The term 5,3 refers to the number of taps in the wavelet fitlers, i.e., 
the number of non-zero (consecutive) values in the basis function support 
for the kernal. Reversible implies that performing the forward and the 
inverse tranformation (using the explicity roimding rules and without 
performing quantization in any form) will result in the exact same number 
at the output as the input. Only a modest and predictable increase in 
precision over the input precision is required for the intermediate 
mathematical terms. That is, there is no systemic distortion introduced by 
mathematical precision. Irreversible implies that very high precision is 
required to ensure no mathematical distortion (exact recor\struction). In 
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practice, however, irreversible filters are almost always combined with 
quatization creating distortion that overwhelms the systemic mathematical 
precision distortion. 

[0142] Figure 24 illustrates one embodiment of a forward transform filter. 
Referring to Figure 24, a highpass filter 2402 is coupled to receive lines 
and the last from line buffer 2401 and generates an output that is coupled 
to one input of the lowpass filter 2404 and is stored in line btiffer 2403. Line 
buffers 2401 and 2403 store one line having a tile width. Lowpass filter 2404 
also receives an output of the highpass filter 2402 from the previous cycle 
that is from line buffer 2403 along with the current x„ line and generates an 
output. The outputs of lowpass filter 2404 for two previous clock cycles are 
delayed through delays 2405 and 2406, providing the filter output one cycle 
in the past and two cycles it the past. 

[0143] Previous outputs of highpass filter 2402 are delayed by delay 2407 
and 2408 such that the current output of highpass filter 2402 and the last two 
outputs of highpass filter 2402 are input to highpass filter 2413. The output 
of highpass filter 2413 is a coefficient in the HH subband and is also input to 
loAvpass filter 2415 along previous output of the highpass 2402 two cycles 
earlier, (two in the past), the output from delay 2408 and the previous 
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output of highpass filter 2413. The output of lowpass filter 2415 is coefficient 
from the HL subbarid. 

[0144] The output of lowpass filter 2404 along with the outputs of delays 
2405 and 2406 are input to highpass filter 2409. The output of highpass filter 
2409 is the LH subband. 

[0145] The output of highpass filter 2409 is also input to one of the inputs of 
lowpass filter 2411 along with the output of delay 2406 and the previous 
output of highpass filter 2409 as delayed through delay 2410. The output of 
lowpass filter 2411 is the LL subband. As the LL subband, the output of 
lowpass filter 2411 is input to a line buffer 2412 the output of which, along 
with the output of lowpass filter 2411, represent the inputs to the next level 
of wavelet transform. The next level of wavelet transform may comprise a 
cascaded version of the wavelet transform in Figure 24. 
[0146] Figure 25A illustrates one embodiment of the lowpass filter, such as 
may be used in transforms (e.g., the 5,3 transform described above) 
described herein. The lowpass filter is designed to create an output based on 
the functions according to the following: 

-x,+2xj-x. 

For the reversible case, the lowpass filter operates according to the following 
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equation: 

X, 

L 2 _ 

[0147] Referring to Figure 25A, an adder 2501 is coupled to add the last 
line with the current line. The least sigruficant bit output represents an 
output of the high pass filter of Figure 25B and is for the irreversible 
transform. The remaiiung bits are input to subtracter 2502 and are 
subtracted from the x^ input to create an output representing the most 
significant bits. These most significant bits are all that is required for the 
reversible case. It should be noted that for an inverse wavelet transform to 
convert filter in Figure 25A into an inverse wavelet filter for use as an odd 
(high pass) filter in an inverse transform, the subtracter 2502 is replaced with 
an adder. Such an example is shown in the high pass filter of Figure 25B. 
[0148] Figure 26A illustrates one embodiment of a highpass filter such as 
may be used in the transforms described herein. For the irreversible case, 
the highpass filter operates according to the following equation: 

4xj-Xq-X2 

[0149] For the reversible case, the highpass filter operates according to the 
following: 
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Xq+X2+2 
4 



[0150] Referring to Figure 26A, adder 2601 adds either the reversible or 
irreversible version of the last line to the current xjine. The output of 
adder 2601 is added to a rovinding term using adder 2603. The rounding 
term may be 2 in the case of reversible or 0 in the case of irreversible and is 
supplied by mux 2602. All except the lower two bits of the output of adder 
2603 is added to the xl line using adder 2604 to produce the reversible 
output. The lower two bits of the output of adder 2603 and the output of 
adder 2604 represent the irreversible output. 

[0151] The use of mux 2602 allows a simple switch to be invoked to switch 
between reversible and irreversible, instead of requiring completely separate 
hardware for both functions or requiring that reversible roimding effect the 
irreversible output. 

[0152] It should be noted that for an inverse wavelet transform to convert 
filter in Figure 26A into an inverse wavelet filter for use as an every (low 
pass) filter in an inverse transform, the adder 2604 is replaced with a 
subtractor. Such an example is shown in the low pass filter of Figure 26B. 
[0153] Figure 27 represents an alternative embodiment of this transform in 
Figure 24 that includes multiplexers (muxes) to perform mirroring at the 
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image boundaries. These muxes include mux 2701-2712. For example, mux 
2701 allows the x^line to be used instead of the last x^ line at a boimdary 
when there is no line in the line buffer 2401 (for example, at the top of a tile). 
Mux 2702 allows the line buffer to provide the other input to lowpass filter 
2404 at times when the bottom of a tile has been reached and there is no 
additional x^ line to be input. Similarly, mux 2703 allows the output of 
highpass filter 2402 to be used as an input to lowpass filter 2404 in cases 
where there is no line in line buffer 2403. Mux 2704 allows the input to 
lowpass filter 2404 to be from line buffer 2403 when there is no output from 
highpass filter 2402. Muxes 2705 and 2706 allow the inputs to highpass filter 
2409 to be the output of delay 2406 and the output of lowpass filter 2404, 
respectively, when an output to lowpass filter 2404 and an output from 
delay 2406, respectively, are not available. The same thing can be said for 
muxes 2709 and 2710, muxes 2707 and 2708, and muxes 2711 and 2712. 
[0154] Figure 28 is a block diagram of one embodiment of an inverse 5,3 
transform. Referring to Figure 28, even filter 2815 is coupled to receive an 
LL coefficient, an HL coefficient, and an HL coefficient from the previous 
cycle from delay 2801. The output of even filter 2815 is coupled to one input 
of even filter 2811, one input of delay 2802, and one input of odd filter 2803. 
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The other inputs of odd filter 2803 are coupled to the HL coefficient from the 
previous cycle via delay 2801 and the output fo even filter 2815 from the 
previous cycle via delay 2802. The output of odd filter 2803 is coupled to 
one input of even filter 2810. 

[0155] A similar arrangement exists with respect to the LH and HH 
coefficients in that filter 2805 is coupled to receive the current HH coefficient 
and LH coefficient along with the HH coefficient from the previous cycle 
output from delay 2804. The output of even filter 2805 is coupled to one 
input of even filter 2811, the input of delay 2806, and one input of odd filter 
2807. The other inputs of odd filter 2807 comprise the HH coeffcient fi:om 
the previous cycle (the output of delay 2804) and the output of even filter 
2805 from the previous cycle (the output of delay 2806). The output of odd 
filter 2807 is input to one input of even filter 2810. 

[0156] Note that the outputs of even filter 2805 and odd filter 2807 are also 
coupled and are stored in inputs of line buffer 2808 and 2809. The size of 
line buffers 2808 and 2809 are equal to the tile width. The outputs of line 
buffer 2808 are input to the other input of even filter 2811 and one input of 
odd filter 2815. The output of line buffer 2809 is coupled to one input of 
even filter 2810 and one input of odd filter 2814. 



70 

[0157] The output of even fUter 2810 is the "C" portion fo the image data that 
is output, is stored in line buffer 2812, and is coupled to one input of odd 
filter 2814. In one embodiment, the size of line buffer 2812 is equal to U the 
tile width. In response to its input, odd filter 2814 generates data 
corresponding to the "A" portion of the image data. 

[0158] The output of even filter 2811 corresponds to the "D" portion of image 
data is input to one input of odd filter 2815 and is stored in line buffer 2813. 
In one embodiment, the size of line buffer 2813 is U the tile width. The 
output of line buffer 2813 is coupled to one input of odd filter 2815. The 
output of odd filter 2815 corresponds to the "B" portion of the image data. 

Other Parallelism Implementation Techniques 

Assignment of Code-blocks to Coders for Parallelism 

[0159] It is useful in hardware implementations to code multiple code-blocks 
in the same tile in parallel. Figure 21 is a memory usage diagram of one 
embodiment of a coder that includes multiple MQ coders, each having an 
associated context model, that may be used to process multiple code-blocks. 
[0160] Referring to Figure 21, each MQ coder is allocated memory (e.g., 
separate memory or some portion of one or multiple memories). In one 
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embodiment, a portion of the allocated memory stores the coded data, with 
the length, zero bitplanes, and coding passes stored in another portion of 
memory. 

[01611 Figures 18-20 show assignment of code-blocks to parallel imits for 
128x128 tiles, 64x64 code-blocks and three transform levels, respectively. 
The assignments are made to balance the amount of coding to be performed 
by each of the parallel coders. In one embodiment, code-blocks are assigned 
so that each MQ coder codes, to the extent possible, approximately the same 
number of coefficients, with a balance between higher level and lower level 
coefficients. Other configurations are possible. 

[0162] Figure 18A-C illustrate embodiments of code-block assignments for 
4:4:4 data when 4, 6, and 8 MQ coders are being used in parallel, 
respectively. In Figure 18C, for 8 xmits in parallel, the code-blocks assigned 
to parallel unit ''H'' (IHH chrominance subbands) will often be heavily 
quantized (have few non-zero bitplanes to code) so it is likely that this unit 
can process more coefficients than the other \mits per unit time. 
[0163] Figures 19A-C illustrate embodiments of code-block assignments for 
4:2:2 data when 4, 6, and 8 MQ coders are being used in parallel, 
respectively. 
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[0164] Figures 20A-C illustrate embodiments of code-block assignments for 
4:1:1 data when 4, 6, and 8 MQ coders are being used in parallel, 
respectively. In Figure 20C for 8 units in parallel, it is expected that units C, 
D and E process more coefficients per unit time than the other imits. 
[0165] The coder of Figure 29 may be used to perform the coding described 
above. For example, each one of N MQ-coders of bit modeling MQ-coders 
2908i.j^ may be assigned to any of A-H shown in Figures 18-20. 
[0166] Note that although even numbers of MQ coders in parallel are 
discussed with respect to Figures 18-20, an odd number of MQ coders in 
parallel may be used. 

Reduced Memory for Storing Coefficients in Hardware 

[0167] Reducing memory usage for storing coefficients when not decoding 
losslessly can make use of zero bitplane information. If the hardware can 
store N bitplanes for each coefficient, decoding can terminate after N 
bitplanes are decoded. Any following bitplanes can be quantized 
(trimcated). 

[0168] Figure 22A illustrates use of a limited nimaber of bitplanes of memory 
for each coefficient during encoding. For example, 8 bitplanes of memory 
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(N=8) can be used to encode coefficients with 16-bits in a normal 
representation. Those coefficients are part of a subband other than the LL 
subband (where the LL subband coefficients are not quantized) that are 
generated as a result of applying a wavelet transform to image data. In one 
embodiment, the wavelet transform comprises a 5,3 wavelet transform as 
described herein. The wavelet transform may comprise multiple 5,3 wavelet 
transforms operating in parallel to generate LL, HH, LH and HL subbands 
in parallel A memory storing coefficients from the wavelet transform may 
be accessed by a context model to perform encoding based on the coefficient 
bits. 

[0169] During encoding, coefficients are stored before the number of zero 
bitplanes is known. Coimters coimt the nimiber of initial zeros for more 
significant bitplanes 8 ... 15. As long as a bitplane 8 ... 15 is all zero, the 
memory stores information (magnitude) for the corresponding bitplane 
0. . .7. Once a one occurs in a bitplane 8. . .15, the corresponding counter 
stops and memory stores information for the corresponding bitplane 8... 15. 
At the end of encoding a code-block, the counters either indicate all zeros for 
a bitplane 8. . . 15 and that the corresponding bitplane 0. . .7 is in memory in 
the case that the coimter stores a value at the end of the memory, or they 
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indicate the starting address for bitplane 8.., 15 data in memory and that the 
corresponding 0,..7 bitplane should be truncated (quantized). Thus, the 
counts act as sideband information to indicate that information stored in the 
memory array from the beginning of a row up to the location in the row at 
the position indicated by the cotint is no longer needed data. The results of 
truncation is that the least significant bitplanes are dropped. 
10170] A separate bitplane of memory can be used to store sign information 
or the sign information can be stored with the significance state. 
[0171] In alternative embodiment, a small amoimt of memory may be used 
for variable length (VL) code information (e.g., run length code) instead of 
counters. This allows a bitplane with a small number of one bits to be stored 
in a small portion of memory for each bitplane. Once the bits have been 
stored in the memory, a context model accesses the memory to obtain the 
bits. However, since each row may potentially contain data that is to be 
quantized, and therefore, need not be accessed and used by the context 
model Figure 22B illustrates a block diagram of one embodiment of control 
logic to control access to the memory. This logic may operate in conjunction 
with or be part of a context model accessing the memory. 
[0172] Referring to Figure 22B, an address, addr, accesses memory array 
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2201 generating a bit. The address and the counter value associated with the 
row of memory containing the address are input to comparison logic 2210. 
If comparison logic 2210 determines that the address is greater than or equal 
to the counter value for the row, then the 1-bit output from memory array 
2201 is output; otherwise, a zero is output* 

[0173] Figure 23 illustrates a portion of memory from a VL code and a 
memory array storing coefficients. The VL code may be used to indicate the 
presence of a 1 bit by indicating an amount of bits to skip until the next one 
is encountered in the row. Thus, this VL code is made to indicate two coimts 
so that the access logic knows where the next bit-plane is. Other VL codes 
could be used to provide more than two cotmts. Use of VL codes typically 
allows one less full bitplane of memory to be used. If the small memory is 
1/32 of the size of the code-block (per bitplane), then a R2[8] code might be 
used. If the small memory is 1/16 of the size of the code-block, then a R2[6] 
or R2[7] code might be used. For more information on R2[8], R2[6]y and 
R2[7] codes, see U.S. Patent No. 5,381,145 entitled ''Method and Apparatus 
for Parallel Decoding and Encoding of Data,'' issued January 10, 1995, 
assigned to the corporate assignee of the present invention. 
[0174] For video where simultaneous operation of the transform and context 
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model/MQ coder is desired, two banks of memory are needed. For still 
image applications, one band of memory is sufficient for sequential 
operation of the transform and context model /MQ-coder. 
[0175] Although the reduced memory technique discussed above is 
described in terms of rows, any arbitrary memory area may be used, such as, 
for example, columns, blocks, pages, regions, etc. Also, separate memories 
may be used. 

Packet Header Processing 

[0176] To create a codestream, such as, for example, a JPEG 2000 codestream 
(or bitstream), packet headers are created. In one embodiment, this 
information may be with a tag tree structure to handle an arbitrary number 
of code-blocks. In certain situations, tile heades for tiles with a limited 
number of code-blocks are created. For example, in the case where a tile 
includes four 128x128 subbands that are each divided into 64x64 code- 
blocks, then there are four code-blocks that are coded together. The packet 
header indicates whether there is any data for a particular code-block, the 
number of zero bit planes if there is data, the length of the coded data and 
how many coding passes the data includes. 
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[0177] Table 9 illustrates one embodiment of a packet structure for packets 
with 2x2 code-blocks and 1 layer. Referring to Table 9, tag trees are only two 
levels high. The references to "z" show where the higher level zero 
bitplanes tag tree information goes, and the locations designated " ^' show 
where the remaining zero bitplanes, coding passes and length infomation 
goes. 



Table 9 - 


Inclusion Information for 2x2 Code-blocks and 1 Layer. 


inclusion 


code 




inclusion 


code 


0000 


0* 




1000 


lllz_000 


0001 


110001z_ 




1001 


lllz_001_ 


0010 


11001z_0 




1010 


lllz_01_0 


0011 


11001z_l_ 




1011 


lllz_01_l_ 


0100 


llOlz.OO 




1100 


lllz_l_00 


0101 


1101z_01_ 




1101 


lllz_l_01_ 


0110 


1101z_l_0 




1110 


lllz_l_l_0 


0111 


1101z_l_l_ 




nil 


lllz_l_l_l_ 



^or 10 or 110000. 

In one embodiment, a 110000 code is used when no code blocks are included 
for ease of implementation. 

[0178] One embodiment of a procedue to write the packet headers for tile 
with a limited number of code-blocks and only a single layer begins with 
initialization that includes the following: 

set minimum zero bit planes, MZP, to maximum value for each 
subband 
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In one embodiment, the maximum value for MZP is OxF for up to 15 

bitplanes or OxlF for up to 31 bitplanes. Larger values may be used for 

implementations that can handle more bitplanes. 

[0179] Next while coding coefficients in each code-block in the packet: 

Save included or not bit 
Save number of zero bitplanes 

If zero bitplanes less than MZP then MZP = zero bitplanes 
Save number of coding passes 
Save length 

The Save included or not bit is set if every coefficient (after quantization) is 
zero, thereby indicating that the code-block is not included. Finally, after 
the information in the tile or subband is processed, the packet header is 
written as follows: 



write "1'' 
for each subband 
write "1" 
first__flag = 1 
for each code-block 

if not included then 
write "0'' 

else 

write "1" 

if first_.flag then 

write MZP in tag tree format 

first_flag = 0 
write zero bitplanes - MZP in tag tree format 
write coding passes 
determine minimum Lblock value 
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write LBlock 
write length 

Note that the LBlock is defined in the JPEG 2000 Standard in section B.10.7.1. 
[01801 Note that the packet header is at least one byte and a JPEG 2000 
compliant decoder can imderstand the information that is written* 
[0181] When there are multiple layers, initialization of the MZP variable can 
be the same as for one layer. While coding each code-block, the included or 
not indication, the number of coding passes and the length are saved for 
each layer. Additionally, the following initialization is preferred. 
first_flag = 1 

initialize Lblock for each code-block 

initialize already included for each code-block to false 

In one embodiment, the LBlock is intialized to 3. The "already included" 

being true means some previous layer had coded data (i.e*, the code-block 

has appeared before). 

[01821 To write the packet header for each layer, the following procedure 

may be used: 

write "1" 

for each subband 

if layer 0 then write "1" 
for each code-block 

if not included then 
write "0'' 

else 
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write ''1" 

if code-block not already included then 
if f irst_flag then 

write MZP in tag tree format 

first_flag = 0 
write zero bitplanes - MZP in tag tree 

format 

set already included 
write coding passes 
determine minimum Lblock value 
write LBlock 
write length 

The "already included" infomation can be a seperate bit for each code-block. 
Otherwise, an otherwise unused value of zero bitplanes can be used to 
indicate "already included". For example, if there are 14 bitplanes, seting 
zero bitplanes to 15 (OxF) can indicate "already included/' 



Reduced Coded Data Not Using ''0" Packets 

[01831 In JPEG 2000, packet headers are roimded up to bytes. In some cases, 
however, a packet header may contain only a single zero bit or a number of 
bits which are less than the number of bits needed to have the packet header 
fall on a byte boimdary. Packet headers are usually roimded to bytes by 
padding. Also, packet header representations are not unique and although 
one typically desires the use of the shortest representation as possible, in one 
embodiment, a representation that is not the shortest possible may be used if 
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the extra bits that are used take the place of bit locations that will have been 
filled by padding. This may be particularly useful in cases where the 
information that is encoded in the extra bits indicates something about the 
next packet in a tile component level partition. 

[0184] For example, if there was a single subband and none of the 2x2 blocks 
was included, one could output a zero packet. However, in the same 
amount of space, one could output a zero to indicate that there is something 
in the packet but there is nothing included in the top level of the tag trees. 
Alternatively, one could indicate that there is something in the tag trees but 
it is 0000 ( or that there is nothing individually in the four). Thus, these extra 
bits may be used to give more tag tree information, which is information 
that would have to occur in the packet header later on and has essentially 
been moved up. By shifting bits up into earlier packet headers, it may be 
possible to reduce the size of the overall codestream by a byte (or more). 
[0185] Whereas many alterations and modifications of the present invention 
will no doubt become apparent to a person of ordinary skill in the art after 
having read the foregoing description, it is to be tmderstood that any 
particular embodiment shown and described by way of illustration is in no 
way intended to be considered limiting. Therefore, references to details of 
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various embodiments are not intended to limit the scope of the claims which 
in themselves recite only those features regarded as essential to the 
invention. 



